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Starting with a covariant spectator quark model developed for the nucleon (N) and the A in the 
physical pion mass region, we extend the predictions of the reaction 77V — > A to the lattice QCD 
regime. The quark model includes S and D waves in the quark-diquark wavefunctions. Within this 
framework it is the D-wave part in the A wavefunction that generates nonzero valence contributions 
for the quadrupole form factors of the transition. Those contributions are however insufficient to 
explain the physical data, since the pion cloud contributions dominate. To separate the two effects 
we apply the model to the lattice regime in a region where the pion cloud effects are negligible, and 
adjust the D-state parameters directly to the lattice data. This process allows us to obtain a better 
determination of the D-state contributions. Finally, by adding a simple parametrization of the pion 
cloud we establish the connection between the experimental data and the lattice data. 
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I. INTRODUCTION 

In recent years, the structure of the baryon resonances 
has been an important topic of both experimental and 
theoretical investigation. Of particular interest is the A 
resonance, the first excited state of the nucleon. New A 
photo-production for large four-momentum transfer have 
been extracted in several laboratories, as Jlab, MAMI, 
LEGS and MIT-Bates 0, i, 0, II H- Simultaneously, 
lattice simulations were performed for both the 7A — > A 
0,0] and the 7A -> A Hi, EH 

transition form factors. 

This last transition can be described in terms of the 
Jones and Scadron multipole form factors [Tl1 | : the mag- 
netic dipole Ml (G* M ), the electric E2 (G* E ), and the 
Coulomb C2 (G* c ) quadrupole form factors. The re- 
action is dominated by the magnetic dipole form fac- 
tor Gy . Although constituen t q uark models (CQM) 
El El EE EE d El, EI H|, with valence quark 
degrees of freedom only, are sufficient to explain some 
properties of the A, they are not sufficient to explain the 
G*m data at the physical point (m, = m% hy = 138 MeV) 
[l2L EH, HH . The inclusion of chiral symmetry and/ or a 
coupling of the quark to the pion field [H, [23|, 0, [2f| 
helps to overcome the limitations of the pure valence 
quark models. The importance of the pion field, or "pion 
cloud" , was also demonstrated using effective field theory 
[H, [27], HH and model reaction mechanisms based on the 
hadronic fields, known as dynamical models [29l [3TI l3l| . 
See Ref. [2l[ for a review of the 7 A — > A transition. 

In prev ious works we successfully described the nu- 
cleon [HEi, the 7 A -> A [34] and the 7 A -> A[HE1 
transitions, by considering a spectator quark model in- 
spired by the Vector Meson Dominance (VMD) mecha- 
nism. When we restrict the A wavefunction to S-waves, 
the spectator valence quark model gives no contributions 



to the quadrupole form factors G* E and G* c [l2j]. Those 
contributions emerge only when D-states are considered, 
consistent with other CQMs [H EH • The D-states im- 
prove the description of the experimental data, but their 
contributions are in general small, with the quadrupole 
form factors at the physical point being dominated by 
the pion cloud [l3[ (see also [H, HH for a review). The 
dominance of the pion cloud at the physical point pre- 
vents an accurate calibration of the valence quark D-state 
contributions and, consequently, the relative amount of 
the contributions from D-waves only and the pion cloud 
is not yet exactly known [IH, [U [3(| . 

The VMD mechanism, which we used to parametrize 
the electromagnetic interaction in terms of the hadronic 
masses, can also be used to extend the covariant spec- 
tator model to the lattice regime. By introducing the 
dependence of the hadronic masses on the lattice pion 
mass we could describe the lattice data for the nucleon 
[HI and for the magnetic dipole form factor G* M [l(| of 
the 7A — > A reaction, with only S-waves in both the 
nucleon and the A wavefunctions (36[. This motivated 
that in the present work we include D-states in the A 
in wavefunction, to show whether the description of the 
lattice data is still possible, also for the G* E and G* c form 
factors. 

In order to include the D-states, and to overcome the 
uncertainty discussed above on their effects, we start by 
realizing that the valence quark effects are expected to 
dominate for large pion masses, exactly when pion cloud 
effects arc expected to be suppressed 



37|. In light of 



this, lattice QCD, in particular in quenched approxima- 
tion with > 400 MeV, where the pion cloud effects 
are negligible (3?j], becomes the ideal laboratory to test 
the valence quark contributions, and to constrain the 
D-states. Importantly, the dependence of the two sub- 
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leading G* E and G* c form factors on the D-states allows 
then the extraction of very relevant information. 

In this work we use the lattice information to de- 
termine more precisely the valence quark effects in the 
77V — > A quadrupole form factors, and consider the 
contributions of the valence quark structure, including 
orbital D-states, in the regime of the lattice calcula- 
tions. First, we started by making a direct application 
of the valence quark model fixed in the physical region 
(tojt = mP hy ) in Ref. [HI], to the quadrupole moments in 
the lattice region. Although the model generates the cor- 
rect order of magnitude of the lattice quadrupole data, 
it fails to reproduce their Q 2 dependence. Therefore, we 
inverted the procedure: we began by readjusting the A 
model parametrization to the quenched lattice data, im- 
posing an overall description of the lattice data for several 
values of the pion mass. From the resulting parametriza- 
tion we generated directly the valence quark contribu- 
tion of the orbital D-states to G* E and G* c at the phys- 
ical point. Finally, by adding a pion cloud contribution 
derived in the large- N c limit, which was established in- 
dependently from our model, we obtain a successful de- 
scription of the experimental data. 

The paper is organized as follows: in Sec. |n]we review 
the formalism associated with the spectator quark model; 
in Sec. IIIII we explain how to generalize the model to the 
lattice QCD regime; in Sec. IIVI we present the results of 
that generalization; in Sec.|V]we present the predictions 
for the quadrupole form factors in the physical region and 
discuss the results; finally in Sec. lVII we draw conclusions. 



II. FORMALISM 

We consider a quark model based on the covariant 
spectator formalism [38l |39| . In this formalism the nu- 
cleon and the A are described as a system of an off-mass- 
shell quark and two non-interacting quarks forming an 
on-mass-shell diquark [l2|, 0, HH, H§] . 

In our model, the quarks are effective degrees of free- 
dom "dressed" by form factors, and the nucleon and A 
wavefunctions are not derived from a dynamical wave 
equation, but given a parametric form consistent with the 
intrinsic symmetries of these systems. By construction 
[T2I |32| the wavefunctions are covariant and reproduce 
their expected non-relativistic limits. The first feature 
is important for the application of the formalism to the 
kinematics of the recent data. 

The nucleon and A wavefunction can be expressed in 
terms of a quark spin (isospin) state together with a spin- 
(isospin-0) diquark state and a spin-1 (isospin- 1) vector 
diquark state [l2 . l32l.l33l ]. multiplied by a relative angular 
momentum state |13j. 

The electromagnetic transition current between the 
nucleon and the A can be expressed [HI, [l3|, [32|, [HI as 



A 



(i) 



where jj is a generic isospin dependent quark current, 
^A; ^iv are respectively the wavefunction of the nu- 
cleon (momentum P_) and A (momentum P+). The 
hadronic current |T]) involves the sum in all intermediate 
diquark polarizations A = 0, ±1, and the invariant inte- 
gral J k = f 2 g rf (27r) a over tne diquark momentum, where 



E s is the diquark on-mass-shell energy E s = y m 2 + k 
(m s is the diquark mass). The factor 3 accounts for the 
flavor symmetry. 

For the quark current we consider the general form 



3 1 =JiY +J2- 



2m 



N 



(2) 



where mjy is the nucleon mass and contain the 

quark form factors, which depend on the transferred four- 
momentum squared Q 2 . They can be decomposed into 
isoscalar (+) and isovector (— ) operators that act in the 
baryon isospin states according to 



h{Q 2 ) = lh + {.Q 2 ) + hh-{Q 2 )rz- 



(3) 



To represent the electromagnetic quark form factors fi± 
(i = 1, 2) we consider a parametrization inspired by 
VMD. In particular, following [H, [H, H2| we used 

h±(Q 2 ) = A + (1 - A) — YT-pri + ca -j^jL (4) 



f 2 ±(Q 2 ) = K± \d ± 



m; 
,2 



i2 1 r>2 



(l-d±) 



(M 2 + Q 2 
Q 2 



M 2 + Q 2 



(5) 



In this parametrization A was adjusted to give the charge 
number density in the deep inelastic limit [32] , m v repre- 
sents a vector meson (m v = m p , m w ), Mh is a mass of an 
effective heavy vector meson simulating the short-range 
structure, and k + (k_) the isoscalar (isovector) quark 
anomalous moments. These two last values were adjusted 
to reproduce the nucleon magnetic moment n+ = 1.639 
and K- — 1.823 [3^]. The coefficients c± and d± were 
adjusted to reproduce the nucleon electromagnetic form- 
factors. In the calculation presented here we took the 
current parametrization corresponding to model II in 
Ref. [13] for the nucleon elastic form factors. The explicit 
values of the VMD coefficients are c+ = 4.16, c_ = 1.16 
and d± = —0.686. For the effective heavy meson we took 
Mh = 2m n. The values for A and the diquark mass are 
1.21 and m s = 0.87 m^ — 817 MeV, respectively. 

For both the nucleon and the A we consider wavefunc- 
tions that contain the correct spin-isospin structure, with 
orbital parts modeled by scalar functions <f> of the quark 
four-momentum squared (P — fc) 2 , expressed in terms of 
the ratio 



Xh = 



{mii — fn s ) 2 — (P — k) 2 
m H m s 



(0) 



where mn represents the nucleon or the A mass (mjv 
and ma)' 



3 



For the nucleon, we consider an S-state wavefunction 
[l2l [32 . 33] , including a mixture of a spin-0 and isospin-0 
with spin-1 and isospin-1 diquark structure. In particular 
we used 



N 



(J) 



where N is a normalization constant. The parameters 
/3i and P2 (P2 > Pi) can be interpreted as Yukawa mass 
parameters. Then Pi accounts for the short range physics 
and Pi for the long range. The corresponding parameters 
can be found in Refs. fl2L fl3L I32I]. 

As for the A, we consider an admixture of S and D 
states, where, as explained in Ref. [li|, the D-state wave- 
function decomposes into a spin 1/2 core (Dl state) and 
a spin 3/2 core (D3 state), both with isospin 3/2: 



N[^ s + a^ D3 + b^! D1 ] 



(8) 



where a and b are admixture coefficients, ^5 represents 
the (symmetric) A S-state, and the remaining two pos- 
sible D states. The normalization constant becomes 
N 



1/vT 



+ b 2 . 



To represent the momentum probability distribution 
of the quark-diquark A system we have, as in [l3j], 



N s 



»D3 



m s {ai 4- xa) 3 
= N D3 
ml(a 2 + Xa) 4 
N D1 f 1 



Am 



2 1(«3 + xa) 4 K + xa) 4 



(9) 
(10) 

(11) 



Similar to the nucleon case, on are momentum range 
parameters, and Nx normalization constants. Since 
with the inclusion of the D-states the S-state can be 
parametrized with only one range parameter (a\), as 
shown in Ref. [l3[, we relabeled all the range parame- 
ters, with ai standing for the average of the values of the 
best model of Ref. In expression (|11[) . the coefficient 
Xdi was chosen to impose the orthogonality between the 
nucleon S-state and the A Dl state (see Ref. [l3[ for de- 
tails). The extra power in Eq. when compared to the 
corresponding equation for the nucleon Eq. ([7]), was in- 
troduced to take into account the G* M falloff observed at 
low and intermediate Q 2 [H, [H, [2l| . The dependence of 
the D-states on the ratio %a was chosen to reproduce the 
behavior of the 7./V — > A transition [l3| in perturbative 
QCD. 

All the wavefunctions are normalized in order to re- 
produce the nucleon fl2l |32| and A pjl, [H, HH charge. 
In particular, one has [l3j: 



/ [<MP,fc)] 2 = / [MP,k)] 2 = i 

Jk Jk 



(12) 



~k 2 4> D3 (P,k) 



k 2 (t> D1 {P,k) =1,(13) 
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Lattice m = 490 MeV 
Lattice m = 41 1 MeV 
Physical point 
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FIG. 1: Quadrupole form factors in quenched lattice QCD 
[lol ]. The lines corresponds to the parametrization of the 
model 4 from Ref. [13 ] . The lattice data corresponds to = 
563 MeV (dashed line), 490 MeV (doted-dashed line) and 411 
MeV (doted line). The solid line is the valence quark contri- 
bution in the physical point (same result as in Ref. [U 



where P represents the momentum in the respective 
baryon frame P = (m^, 0,0,0). The variable k defined 
as k = k — ^r-P was introduced in Ref. [l3[ to represent 
the D-states. 

In summary, the nucleon (S-state) wavefunction con- 
tains two range parameters (/3i and fa), the A S-state 
wavefunction a single parameter (ai), and the A D-states 
contain three (a,, with i = 2,3,4) parameters. The D- 
state admixture is controlled additionally by the a and b 
parameters defined in Eq. ([5]) . All the parameters associ- 
ated to the quark current and the nucleon wavefunction 
were fixed by the nucleon elastic data at the physical 
point in a previous work (32j. They are kept unchanged 
here, where we fit only five parameters to the lattice QCD 
data: the three range parameters and the two admixture 
coefficients for the two D-states. 



III. EXTENSION OF THE SPECTATOR MODEL 
TO THE LATTICE REGIME 

Now we consider the extension of the model presented 
in the previous section, and based on the valence quark 
degrees of freedom, to the lattice QCD data region. As 
the pion cloud effects are expected to be suppressed for 
TOtt > 400 MeV [I?]], it is justified in this region to con- 
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JTItt TTtjV 771A Tip 



0.563 


1.267 


1.470 


0.898 


0.490 


1.190 


1.425 


0.835 


411 


1.109 


1.382 


0.848 


0.138 


0.939 


1.232 


0.776 



TABLE I: Masses in GeV, considered in lattice QCD [T(| and 
at the physical point. 



sider the valence contributions only. Therefore, we gen- 
eralize the original (valence) spectator model by consid- 
ering the implicit dependence of all the mass parameters 
on the pion mass. Since the wavehmctions depend on the 
ratio defined in Eq. ([6]) the diquark mass scales out from 
the current and consequently from the form factors [321 ] . 
The dependence of the model on the pion mass appears 
through both the quark current and the baryon wave- 
functions. The quark form factors (Jll)-© in the current 
operator (fSJ) are written in terms of a VMD parametriza- 
tion, and therefore depend on a vector meson mass m v 
and an effective heavy meson mass — 2m jy. The 
extension to lattice is done by considering the nucleon 
and p masses from the lattice calculation as an implicit 
function of m^. 

The nucleon and A wavefunctions ([7]) and (|5|- ([TTj) are 
represented in terms of the (adimensional) momentum 
range parameters and o^) and the kinematic ratio 
Xh of Eq. ([6]) . Because the mass dependence of the wave- 
function enters in this ratio, we expect only a weak de- 
pendence on the range parameters (cti and Pi) near the 
physical region. In fact, this weak sensitivity of the range 
parameters to the pion mass was already verified in the 
work of Ref. [36[ for light pions. For this reason we use 
the same range parameters both in the lattice data re- 
gion and at the physical point, neglecting any pion mass 
dependence of the range parameters in the considered re- 
gion (m^ < 600 MeV) . For heavier pions, the interaction 
becomes almost pointlike, and at least the short-range co- 
efficients should be corrected, or an explicit dependence 
on the pion mass introduced. 

The procedure presented here corresponds to the case 
of Ref. [3y| with M x = +oo (M x is the constituent quark 
mass in the chiral limit). A more detailed treatment is 
possible with finite values of M x , but for the quadrupole 
lattice data the corrections are small when compared 
to the statistical errorbands and the differences between 
the datasets associated with different values for the pion 
masses. 

To start with we took the best valence quark model 
presented in Ref. [13|, fixed for the physical data case 
(tott = mP hv ). The analytical expressions for G* M , G* E 
and Gq are presented in Ref. [13|]. No readjustment of 
the parameters of the quark current and of the nucleon 
and A wavefunctions was made, except for the nucleon, A 
and p masses, which are explicit functions of m T , as in the 
lattice calculations. The mass parameters are presented 



m n (GeV) X 2 (Gl 4 ) 


X 2 (G* E ) 


x\G* c ) 


x 2 


0.563 


0.569 


0.483 


0.853 


0.618 


0.490 


0.544 


0.290 


0.668 


0.513 


0.411 


1.956 


0.548 


1.163 


1.406 


Total 








0.842 



TABLE II: Quality of the quenched fit in \ 2 (partial and 
total) for the three form factors at the respective pion mass. 
Quenched lattice QCD data from Ref. [Tol ] . 



a 2 , «3, «4 a, b 
0.3421 0.3507 0.0856 
0.3377 0.0857 

TABLE III: D-state parameters of the A wavefunction as re- 
sult of the fit to the quenched lattice data 10]. The coefficient 
Adi = 1.0319, in Eq. (|ll|l is determined by the values of 0:3 
and Q4. 

in Table [J 

Figure [T] shows the results of model 4 of Ref. [l3[ 
extended to the quenched lattice QCD calculations of 
Ref. We only show G* E and G* c , since G* M does not 
change much from 36]. The conclusion is that we can- 
not reproduce the lattice data accurately, but the pre- 
dictions of the model have the right order of magnitude. 
The poor description of G* E and G* Cl obtained from tak- 
ing the physical model to the lattice domain, contrasts 
with what happens for the dominant G* M form factor ob- 
tained in Ref. [36| . It shows that the lattice data is more 
sensitive to the D-wave components of the A wavefunc- 
tion, crucial for the quadrupole form factors, than to the 
S-wave components, that dominate G* M . This makes the 
lattice data for the quadrupole form factors extremely 
interesting, as a potential source of information on the 
D-wave effects in the hadronic structure, and thus indi- 
rectly, on the magnitude of the pion cloud effects. In 
Ref. [ijl, when the model was applied to the physical 
point, the valence quark contributions were not explicitly 
separated from the pion cloud contributions. For that 
reason the estimate of the valence contribution could not 
be made very precise, as we confirm in this extension to 
the lattice regime. 



IV. ADJUSTMENT OF THE D-STATE 
PARAMETERS TO THE LATTICE DATA 

Given the results of the last section, we decided to 
change our strategy in the process of comparing the quark 
model results with the lattice data. Instead of trying to 
use the quadrupole data in the physical region, where va- 
lence and pion cloud contributions are both important, 
we fit first the lattice data and extrapolate to the physical 
point. With this procedure we avoid ambiguities related 
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FIG. 2: Magnetic dipole form factor G* M in quenched QCD [Io| . The valence quark contribution at the physical point (solid 
line) is included as reference. 
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FIG. 3: Best fit of the quadrupole form factors in quenched 
lattice QCD [l(J. The lines have the same meaning of the 

Fig.m 



to the exact contribution of the pion cloud mechanisms 
and their entanglement with the D-wave effects. The pa- 
rameters which are varied in the fit are the admixture 
coefficients a and b and the three range parameters on 
(i = 2,3,4), all associated to the D-state scalar wave- 
functions. We fit these five parameters to the lattice 
data. The S-state parameter (a± ) was not fitted, and the 
value ol\ — 0.3366 from Ref. [l3| was used. 

The partial and total results for the x 2 obtained are 
presented in the Table [ill The final results for the ob- 
servables are presented in the Figs. and [3] The D-wave 
model parameters associated with the fit are presented 
in Table [ml 

The obtained range parameters are larger than the 



ones from the fit to the physical data [TJ], suggesting 
that the D-states are less peripheral (i.e., have a shorter 
range in configuration space) than inferred in Ref. 
where we used an indirect estimate of the valence quark 
contribution at the physical pion mass point, based on 
specific assumptions about the pion cloud contribution. 
Another interesting point is that the D-state range pa- 
rameters cti — 0.338 — 0.351, with i — 2, ..,4, do not 
spread over a large region, suggesting one single value, 
«i ~ 0.344, as the signature range of the D-state regime, 
slightly larger than that of the S-state range ot\ ~ 0.337. 
We note that, in this case, it is the additional power in 
the D-state scalar wavefunctions (see Eqs. (|9l)- ([TT]) ) that 
implies a more peripheral character for those states, when 
compared with the S states. Indeed, as for low k, x = 

[39l | . a higher power in momentum space corresponds to 
a more peripheral effect in configuration space, as it is 
to be expected from a D-wave contribution. The best fit 
corresponds to an admixture of 0.72% for both D3 and 
Dl states in the A wavefunction. 

Compared to the results obtained by using the physi- 
cal data alone, the significant change occurs for the per- 
centage of the Dl component, that drops from 4.36% to 
0.72%. As for the state D3 the differences are minor. In 
Ref. [l3[ the percentage was 0.88%. As a ~ b, we con- 
clude that the quenched lattice data is consistent with an 
equal admixture for both D-states. The initial number of 
effective parameters needed for a good fit can be reduced 
from five to four. 

In Fig. [2] the quenched G* M data is very well described 
by our model [which is also stated in the Table |TT] in 
the column x 2 (G*m)]- The exception is the case for the 
lightest pion mass m, = 411 MeV, where pion cloud 
effects may start to be important [37] but are absent in 
the valence quark model. 

The model describes fairly well the lattice data for the 
quadrupole form factors. The quality decreases for the 
lightest pion mass, which is due to the omission of explicit 
pion cloud effects in our approach. It is encouraging that 
the x 2 values are lower than the ones found for the fit 
in the physical region [l3j , indicating that the procedure 
used here is more natural. Still, it should be said that 
the % 2 's obtained also possibly reflect the still poorer 
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quadrupole lattice statistics and the narrower range of 
the lattice data, when compared to the experimental data 
or even to the lattice data for G* M . 

A study of the dependence of the jN — > A form 
factors on the pion mass and Q 2 was also considered 
The lattice QCD data for G* M in 



in Refs. [2€ 

Refs. 0, [26[, manifests a significant difference from the 
more recent analysis of Ref. 10]. For similar pion mass 
the results of the Ref. are larger than the ones pre- 
sented in Ref. 11011 - 



V. QUADRUPOLE FORM FACTORS AT THE 
PHYSICAL POINT 



After the parametrization of the D-states was obtained 
from the lattice data, we can now apply it to the physical 
region. This requires that we use the experimental data, 
in the form of the two ratios 



Rem 



G* E (Q 2 ) 
g *m(Q 2 )' 



R 



EM 



|q| g* c (Q 2 ) 

2m A G* M {Q 2 ) 



(14) 




where |q| is the photon momentum in the A rest frame, 
and we use the empirical parametrization of Ref. [271 ] : 



G* M (Q 2 ) =3G D exp(-0.21Q 



In the last expression Gd = (l + Q 2 /0.7l) represents 
the dipole form factor. The quality of the parametriza- 
tion was studied in the Ref. [131 ]. 

Since D-states extracted from the lattice data and ap- 
plied to the physical region include only the contribu- 
tion of the valence quarks, they necessarily underesti- 
mate the experimental data. To fill the gap between the 
valence contribution and the experimental data we have 
to consider contributions from the pion cloud. In partic- 
ular we consider the pion cloud parametrization used in 
Ref. [131 ] , where the pion cloud contributions to G% and 
G* c were determined using large- N c relations E| be- 
tween those form factors and Geu (the neutron electric 
form factor): 



GUQ 2 ) 



m N \ 3/2 m 2 A - m 2 N G En (Q 2 ) 



GUQ 2 )= l2mN 



mNinA- 



2V2 Q 2 

Geu{Q ) 



Q 2 



(16) 
(17) 



To evaluate G_e„, we took model II of [32j for the nu- 
cleon. The results are presented in Fig. [4j In that figure 
we compare the final results for G* E and G* c and include 
the lattice data to show the magnitude of the valence con- 
tributions. The valence contributions are also compared 
with the parametrization of the valence contribution from 
the Sato and Lee model [3l| . Note that the Sato and Lee 
parametrization gives a contribution similar to our model 



for Q 2 > 0.5 GcV 2 , for both G% and G* c . It lies above 
our results, overpredicting the lattice data, for lower Q 2 . 

In conclusion, by fixing the D-state components by the 
lattice data and considering a pion cloud parametriza- 
tion, derived from the large-A^ limit at the physical 
point, we obtained a fairly good description of the 
quadrupole lattice data, in the range Q 2 < 1.5 GeV 2 . 
The exception is the region Q 2 < 0.2 GeV 2 where a 
small Dl-mixture, when compared with Ref. [LSI ], under- 
predicts the G* c data. Note however that there is some 
discrepancy between different experimental data in that 
region [13]. The planned data from the CLAS collabora- 
tion for that range would be important to clarify the low 
Q 2 behavior of G* c [SI S| . 

A complete lattice QCD dataset is available in 
Ref. [l(J. It presents lattice data in the quenched ap- 
proximation, and the unquenched data based on Wilson 
and also on a hybrid action. In this work we restrict our 
application to the quenched data. There are three main 
reasons for this restriction: 

• There is a significant discrepancy between the 
quenched and unquenched data, particularly for the 
results of G* M with heavier masses. In this regime 
we would expect small pion cloud effects, imply- 
ing negligible differences between quenched and un- 
quenched results. In Fig. [5] we compare the lattice 
data corresponding to = 563 MeV (quenched) 
and ra„- = 594 MeV (hybrid action). There is a 
significant difference between those two data sets, 
with the Wilson data associated with m n =691 
MeV being more consistent with the hybrid action 
(m n = 563 MeV). 

• There arc differences between the two unquenched 
results, in particular for the value of m p . The ex- 
tension of our model depends on the (quenched) 
p mass. It is not clear whether the extension of 
our model is justified for the unquenched calcula- 
tions, where the nucleon, A, and p masses would 
differ from the quenched masses. We would expect 
only minor differences for heavier pion masses (say 
m w > 480 MeV). However, the significant difference 
between the p mass for the Wilson action data with 

= 509 MeV (m p = 887 MeV) and the hybrid 
action with = 490 MeV (m p = 949 MeV) is 
difficult to explain. 

• Finally, there is only a limited number of un- 
quenched quadrupole data points for large pion 
masses. Since the unquenched lattice data for light 
pion masses as 353 MeV (hybrid action) and 384 
MeV (Wilson action) are expected to be contami- 
nated with pion cloud effects, which cannot be sim- 
ulated by our valence quark model, those points 
would have to be excluded. We would then be left 
with 6 or 8-9 quadrupole points, respectively, for 
the Wilson and hybrid action (to be compared with 
21 from quenched data), and with such a small 
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FIG. 4: Extension to the physical region using the quenched parametrization. The dashed line represents the valence contri- 
bution. The solid line represents the combination of valence and pion cloud effects. Physical data from Jlab P, Q|, MAMI 
LEGS [3] and MIT-Bates Q. Quenched lattice QCD data from [Tj}. Sato and Lee parametrization from [3ll |. 



number of constraints the fit would naturally be- 
come meaningless. 



Apart from the disagreement observed in the descrip- 
tion of the magnetic dipole form factor, the Wilson and 
hybrid action lattice data suggest a weaker falloff of the 
electric quadrupole form factor G* E when compared to 
the quenched prediction. This result is also observed for 
the physical point extrapolation. 

Once the differences between the two unquenched re- 
sults are understood, and the disagreement between 
quenched and unquenched results for ~ 600 MeV is 
clarified, it would be interesting to use also unquenched 
data to extract the contribution of the D-states, using 
the procedure suggested here. The increasing of statis- 
tics in both quenched and unquenched lattice data could 
also help to constrain the effects of the A D-states in the 
jN — ► A transition. 



VI. CONCLUSIONS 

In this work we study the valence quark contributions 
to the 7^ — > A transition in the lattice QCD regime, 
in the framework of the covariant spectator formalism. 
The nucleon and the A wavefunctions are not derived 
from a wave equation but are parametrized in terms of 
the nucleon and A symmetry structure for spin, isospin 
and angular momentum. By construction, the formalism 
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0.5- 



m = 563 MeV (Quenched) 
♦ m* = 594 MeV (Hybrid) 
■ m =691 MeV (Wilson) 



•t • • . s. 



0.5 



Q 2 (GeV 2 ) 



1.5 



FIG. 5: Dependence of G* M with m,n, for quenched (m n = 
563 MeV), Wilson (ra, = 691 MeV) and Hybrid (m, = 594 
MeV). The lines correspond to the result of the VMD model 
with mjv, m p and ijia associated to the quenched data for 
= 563 MeV (dashed), Wilson action: = 594 MeV 
(doted) and Hybrid action: = 691 MeV (dashed-dotted). 



includes only valence quark degrees of freedom, whereas 
meson cloud mechanisms are not taken into account. 

As discussed extensively in the literature (see Refs. [l3l . 
|26|). valence contributions do not dominate the 
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quadrupole form factors at the physical point, where the 
pion cloud effects dominate instead, while the opposite 
happens in the lattice QCD regime. In an attempt to ex- 
plain the significant difference between the experimental 
data and the emerging simulations of lattice QCD with 
decreasing pion mass, but still small chiral effects asso- 
ciated with the light pions, we started by comparing our 
quark model directly to the lattice data. In contrast to 
the dominant contribution controlled by the nucleon and 
A S-states, the A D-states, being a second order correc- 
tion, are more sensitive to the lattice data. This sensitiv- 
ity provides a clean evaluation of the valence quark con- 
tribution, since in the physical region the valence quark 
contribution is masked by the overwhelmingly larger con- 
tribution of the pion cloud, and consequently the D-state 
paramctrization cannot be accurately constrained. 

Accordingly, we found that by fixing the D-states by 
the physical data first does not lead to a good descrip- 
tion of the lattice QCD data. We also verified that, in- 
versely, when the D-states are first fixed in the lattice 
QCD regime, then a good description of the physical data 
is possible. In fact, adding the valence quark contribution 
extrapolated from quenched lattice QCD to the physical 
mass regime, with an estimate of the pion cloud based 
on the large- A^ c limit [HI, H(| , the experimental data for 
"fN — > A quadrupole form factor is well described. An 
even more accurate description of the experimental data 



can in principle be obtained by considering more and 
more precise lattice QCD data, and a more sophisticated 
estimate of the pion cloud. 

The fit to the lattice QCD data varied five parameters 
associated with the valence D-state states. The result of 
the fit suggests an identical admixture of the Dl and D3 
states. 

We conclude that lattice QCD data are important to 
constrain valence quark models. Since for lattice calcula- 
tions with m T > 400 MeV valence quark effects dominate 
over the pion cloud, lattice data can be used to study and 
separate those effects. 
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